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[Jan. 6, 1836.] 
[1. The principal function for an elliptic orbit.] 


For the undisturbed motion of a planet, with Newton's law of attraction, let* 


$ (£2 i2. . 
0 2 7; 


£i, 7 Ti being the two rectangular heliocentric coordinates and the heliocentric distance of а 
planet, and u being the sum of the two masses; so that the differential equations of motion of 
the second order are 


Я = — шг, qur? 
and, by my general theorem of dynamics,t 


бе = (3, 686+ mom (E- 827 ar 


in which the coefficient of ôt is constant, 
E PERHAP) = E p E+). 
t To 


The function s will therefore satisfy the two partial difterential equations following:t 


Lo e Cae И 
5258) "2\8 п’ 86 2V96)  2Vom EN 


в being considered as dependent on the 5 independent variables £,, ту, £,, т; and t, besides the 
constant p. 


* [See Appendix, Note 5, p. 624.] 
T [Cf. pp. 160, 166.] 
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Let r,2«-- B, rg2«— B, (£j —£9)* - (q;— 99)? = (r--rj)?sinA?. I have found,* first, that s 
depends only on the three variables «, À, t, being independent of B; and, secondly, that when А 
is given, 8/У шх depends only on t V/u/a?; so that if we put 


um EOM 
в=к\у ux, ser 


к will be a certain function, of which the form is to be discovered, of the two variables А and v. 
And if the motion is nearly circular, these two variables А, v will be nearly equal to each other. 
(In some loose sheets of investigations begun in May, 1834, I have used the symbol v to denote 


2 


Taking the variation of the expression s — x V ра, we find 


The -Vua (ОВА ае E ЧЕГ. ба; 


in which 
A _ £j0£,-- 9n, 686+ тед 
Vl. Licini dryer o Apre cp | 
бу=— X WC (icc, f Eo Éo + Лобо 
5 2r, 2rg ч 
апа 
(£, — Eo) (8,— 60) + (7: — то) (бт, — 999) — (r + ro) sin А? (òr Tn) 
ФА = O E о 1 э м = 
(r + rg? sinAcos A 
Hence 
_ (к [ш Зв бк _ 24 (набат. Eo Eo + подто 
yd ($ J! 4 о? ду илә 8A 2r, 2% 
совес 24 [ш дк P A ay 
Ыы 9 p" P qua Eo) (8&— 80) + (m — 99) (бт — 559) * 5 t; 
therefore " 
à 8 
SA Bi) у =АЛ+Воу—лЫ), 
AB "^ B (&— £o), ире В (qi — то), B C, 
if we put, for abridgment, 
4 "(t аа А8к Зудк ЈА, pAs [ы (у px. 
-(- 2 ÀA 48» 2 8а ЗАМа’ 2a v! 
and since 


E-e GT- чье 


To 
7} — 15 + 4o? sin А? 
т; : 
rà — ri + 4o? sin А? 


- но Po у. „4... д 


To 


7 (£ — £y) + m (99) = 


* [See Appendix, Note 5, p. 624.] 
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the two partial differential equations become, when multiplied by 7;, rọ, and then added and 
subtracted, and divided by 2«, 28, 


E= C 4- 1A? + 2B?o? sin А + 24 Basin X2, 


0=С+142 + 2 В?о? sin Y? + 24 Вх. 
Resolving these two equations for the two quantities C + $4? + 2 В?о? sin А? and 24 Bx, we find 


a? = 2 sec А?; 


: (2C + A? + 4B?a? sin А?) = E 
therefore 

(4 + 2Bssin Xy 2 [вом C — "sin Ано 3°), 

(A —2B«sin А) = 2 (8м С+іпдввод?); 


and, substituting for А, B, C their values, the equations become 


(«2л 0-3 x) e | sec À? + 162* — 32sec X, 


5А "^8 8А 


бк g „Эк дк 
Bus у 797 


(«- 3,2 8). Е) +16 (8-2) 0.] 
These two equations ought to be compatible; that is, they ought to give equal values of 
2 
by elimination y and differentiate it with respect to А. Without this previous elimination, we 
may differentiate the two partial differential equations separately and successively with respect 
to Лала v, and then eliminate, between the four resulting equations, the three partial differential 


ipia 
coefficients of the second order e es Жубу? т 4 the resulting equation ought to be true either 


(k-2tana gi 3 сао. 


[Непсе 


whether we calculate by elimination x and differentiate it with respect to v, or calculate 


identically or at least in virtue of the two partial differential equations of the first order. 
After some reduction it is easy to see that this condition of compatibility is 


бк [дк бк\? бк 
— 2 а verg ВС 2 2 
0= – 45, >) + 16 tan А? (4 sec à? — 9+2 (5) в Ча sec А (>) +4 
\ ёк |? бк\? 
e 2 е Om 2 
Sec А (5) (sec n(x) 8sec À |, 
in which, by the first partial differential equation, 


бк [дк 4 òr \4 " бк 2 а 
-45 (5). = 16 tan А? + 5) sec À -8 (2) sec А?, 


во that the condition is in fact gatisfied, and the two partial differential equations are compatible. 


28-2 


WWwW.rcin.org.pl 


220 V. NEARLY CIRCULAR ORBITS [2 


[2. Approximations to the principal function in the case of nearly circular orbits.] 


Of these two partial differential equations, the second, namely 


бк дк\ бк 


seems well adapted for determining the coefficients of a series such as the following: 
(v — А)? (v—A 


K — Kg Ki (v — À) *- kg а 779689 


lu 
+ etoe., 


if it is possible, as it seems likely that it is, to represent к by such a series for the case of nearly 
circular orbits.* Assuming this series, we have 


8 , , – А)? , — А) р 

к =кр— ку (P-A) (4) А (i) EE (1) ete 
ò —A9 — \)З 

(А) + ә ) s, eMe к, + ete.; 


к; being a function of А only, and к; being its first derivative with respect to А. Thus, comparing 
the terms independent of v — À, we find the equation 

0 — («y — 2 tan А (ко — к!) — ЗАк!} (ко — ку) + 8tan А; 
in which, from the theory of circular motion, ко = ЗА, and therefore к, = 3, so that the equation 


becomes 
0 — (к, – 3) {к, (3A — 2tan à) + 3 (2tan A — 4)) + 8tan A, 
that is, 
0 = кз (ЗА — 2tan à) + 12«, (tan A — А) + 9А — 10 tan À 


= (i, — 1) {к (ЗА — 2tan А) — (9A — 10 tan 2); 


therefore either к, = 1, or else 
9А— 10%апА. 


— 3À-2tanÀ' 
and in order to decide which of these two roots we are to take, it seems necessary, in the present 
method, to employ the other partial differential equation, namely, 

дк дк 


абе аер ы үө И 
-(«- n A Keg + SA sec 


This equation gives, by making у= А, 
0 — {ко — ЗАк, + 2 tan À (ку —к)} + 4 (к, — ко) sec А? + 16 (ie, — 2 sec А); 
therefore, since ку = ЗА, ку = 3, 
0 = {(3А — 2tan А) к, + 3 (2 tan A — 4))? + 4 (к, — 3)? sec А? + 16 (к, — 2вес À?) 
= к? {(ЗА — 2 tan Л)? + 4 sec А) + 2«, (3 (ЗА — 2 tan À) (2 tan À — A) — 12 sec А? + 8} 
+9 (2 tan À — A)? + 4sec À? 
= (к, — 1) {к, (93? — 12A tan À + 8 tan А? + 4) — (93? — 36A tan À + 40 tan А? + 4)); 
* [r.—»—2 is small when the orbit is nearly circular. See Appendix, Note 5, p. 626.] 


Ki 


8к 


шай 2 
55 эвес 3°). 


WWwW.rcin.org.pl 


2] V. NEARLY CIRCULAR ORBITS 221 
therefore by this equation we have either к; = 1, or 


9 (2 tan À — Л)? + 4 sec А? 
К. (31 — 2tan А)? 4вес À?’ 
9A— 10 tan А 


and, since, by the former equation, we had either к, = 1 or к, = '34-2tanÀ ' 


we must conclude 
that к, =1. If we then neglect (v — А)?, we have simply 


| к= Kg-F k (УА) = 2 - v, | 


Since the two partial differential equations separately gave, as we just now said, ambiguous 
results, but such that the ambiguity was removed when the results of the one equation were 
compared with those of the other, it seems evident that some combination of the two equations 
must be capable of giving unambiguous results, at least for the coefficient к, , and probably for 


the other coefficients. Accordingly, if we put S А under the form D "u е е , and then eliminate 


because ку — ЗА, к; = 1. 


the square of that г which is not added to ar , we shall accomplish the purpose proposed; and 


the same object may be accomplished still more simply by introducing a new variable в, which 


shall be equal to v — А, and shall therefore be small for nearly circular orbits, and thus changing 
бк, ёк бк, дк дк ak? 
5 фо — P and =— SA to = 3A x in the two partial differential equations, and then eliminating E 5, 1 


between them. 


In this manner we have 


E T 
M ad? , + жа , 13 + ete 
о "Mur Мыз a 
бк 2 e 


к; being some function of А, and к; = x ; also the two partial differential equations become 


дк бк дк\? А бк 
0= fk- 2tana èr x (2tanA- maea) ala) sec À +16 (S 252028), 


бк) [к дк 
0=|к- 2anAP + (2tanà— за)“ EIS) tana 


2 
The coefficient of (>) in the first is (2tanA— ЗА — 34)? + 4sec à?, and in the second it is 


2 
— (2tan A — ЗА — 3:); and since the terms of the form const. x « (>) cannot cause ambiguity, it 
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is sufficient for our present purpose to multiply the first equation by 2tan A— ЗА and the second 
by (2tanA— ЗА)? + 4sec 5?, and to add the products together. In this manner we get 


бк бк бк бк бк 
0- (к 2tanAg + (2anA- 83) 3 | (2tanÀ- 3) [e 335 —® | 


бк дк А бк дк бк k 
+4 (5-х) sen к-за |+ 16 (S - 260A | (2tan A — 3) 


+ 8tan À {(2 tan À — ЗА)? + 4 sec А); 


that is, 
дк бк à бк бк дк 
o- [(«- (2tan A — 32) 4-2 (2-- 33 tan А) с + (9) — I2ManA-4)s] (ni -*) 


+ 16 (2tan A — ЗА) (Ss — 256029) + 8tan A((2tanA — А+ 46024); 


бк дк\? дк бк\ [ёк дк 
o- (.- Ag - а) (2tanA- 3) В («9A — 96) (8-9) 


+ 16 (2tanA— 33) + 8tan À (9А — 4) + 96А. 


tazen. = 0, we get 
(А 


Making, in this, к= ко, A ed, E 
0 = (ко — 35)? (2 tan A — ЗА) + (9A? — 12A tan À — 4) (xo — ЗАко) (ку — ко) 
+ 16 (24ап А— ЗА) к, + 8 tan À (942 4) + 962; 
and if we further employ the values ку = ЗА, к= 3, we find 
0 — 362? (2 tan À — ЗА) — 6А (933 — 122 tan A — 4) (xı — 3) 
+ 16 (2tan А— ЗА) к, + 8 tan A (9? — 4) + 96A 
= 2 (ку — 1) (92? -- 4) {4 tan À — 3А); 
therefore к, = 1, as before. 


In general, 
дк дк)? 92? + 4 дк x) (- бк >) 

o-[- eem) (аи ад] а И” га 
бк А 93-44 | 
+16 5; +4 (9А —4)4 122 ау) 

дк дк 

а formula which, when „= 0, к= ЗА, к 3, à" becomes | 
91? 4 "d 99-4 
0 = 36A? + 6А (o soia) (к, — 3) + 16, + 4 (93? — 4) + 12A Ly wap id 


and gives x, = 1. When we make к = ЗА + 0, = =3, «ode 1 + кое, and neglect :?, the formula becomes 
L 
932+ 4 


214 А 
O= (2) (eie EE s) С? + ки) (92929 6 161-9) 4 (94 4) Age gy 


that is, бй 
0= 24A + (62uc — 4) (4+ i) + 16к, 

9%? + 4 99-4 
= к [603+ 6 стт АЗА тта; 
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that is; finally; — ини 
Thus 
12 
в, 
neglecting (3. 


If we make the same substitutions, («= ЗА +, E 8, - 1-4 Ks ) , and still neglect 12, in 


the partial differential ira y 


бк дк 
0 = fk- Зап x (2tanA- 3A— 395. (5% Sc) + вал, 
we find 
0 — (3A 4- .— 6 tan À + (2tan А — ЗА — 30) (1 + «3:)) (2 — kat) + 8tan А, 
that is, 
0—2(1— 3-4 (2tanA— 3А) к,) + 4«, tan, 
or 
E Млн вя 
Ка 41tanÀ—3À 


аз before. But it is remarkable that this value has been deduced without ambiguity from the 
tolerably simple partial differential equation which we set aside before as likely to give only 
ambiguous results. 

If we return to the system of variables к, А, v, ibtd of к, А, в, and therefore resume the old 
partial differential equation 


бк „Эк бк 
0 = (-- 2tanA, 3 "s 5. + Stan 
together with the old development 
: (и А)? P ых 
K — Kg tk (v — À) 4- ко Tk Ks + eto., 


2 


we may propose to determine the several coefficients ку, кү, кз, ... as the values of the partial 
differential coefficients 
oue rue gw 


бу’ à, бы’ 
when v becomes equal to A. Thus, differentiating the equation with respect to v once, we find 


o2 - BO ава аа в кь 


E SÀ п o ih УЗ. ma A ETEM 
Hi бк 8°к бк БЫЗ ёк ёк _ 52к мы Atana Se 9*« 92к 
"ЧАБА ГУ А) | 5ABv 55, 555» 
Sk бк iM: CS E 


which, if we make у= А, к= ко, =K, becomes 


БАС КОТ Е то О УУ бао о gA 
0= ЗА{(ко– ку) ка + к; (1 — кз)} + Эк; (ко — K1) — ко (кү — Ka) + 4 tan А (Ko — ку) (ку — кз), 
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that is, 0 = ((3A— 4 tan А) («o — к!) — (ЗАк, — ко} (kg — ку) + (3Ак + 21) (s — 1); 
and this, when we make ку = ЗА, ку — 1, kọ = 3, к = 0, becomes 
Ka = sie 
?" 4tanA— ЗА’ 
as before. 


Accurately, for elliptic orbits, the relation between к, А, « may be found by eliminating the 
auxiliary variables e and v between the three following equations:* 


3v +esinv { 


2? А Пед» v—esinv 
jp ES (1 —ecos v)!" 


АЕ 
(1—eoos v)? 


and if the orbit be nearly circular, then г and e will be small and of one common order. 


Moreover, if a be the semiaxis major of the elliptic orbit and « (as before) the semi-sum of 


the two given distances r and rg, then 
a дк дк 


a 5 5’ 
А being treated as a constant. Hence, in series, 
х | 12 13 
ака + at Ks ао atete; 


and hence we might have easily foreseen that к; must be equal to 1, because it is the value of 
оја for circular motion. The finding of this value к; = 1 was almost the only difficulty in the use 
of the partial differential equation 


бк дк\ бк 
o- (к-з "s 3.) gs * Stand; 
and if we assume, as we now see that we may, 
—А y—AÀ 
k22)4 v4 ica га АУ еке, 


the labour of calculating кз, кз will be but trifling and the operation will become one of elemen- 
tary facility. In fact, to resume, we have, neglecting 13, 


бк А 12 
5=2-— ка + (ка — кз), 
бк e 
м ка, 


° 
к=ЗА+ +, и=А+ к 


and, if we substitute these values in the partial differential equation, we find 


к 2tanA? За = — 4tan À + ı {r (2tan A — ЗА) — 2} 


2 
+ á {кз (2tan A — ЗА) — 5x, — 2x5 tan А). 
* [Cf. Appendix, Note 5, p. 625.] 
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2 
Therefore, multiplying this by = =2—к„— 5 (кз— 5), adding 8 tan А, and equating to zero the 
coefficients of ı and ‹?, we have 
0= к, (4$ап^— 3А) – 2, 0= к (4%апА– ЗА) — кз (2tan А— ЗА) — Зк — 4r; tan А; 
the first equation gives 
2 
Ка 4tanÀA— 3À' 
as before, and therefore 
‚_ —2(1+4%ап А?) 
2 (4tanA- 30)? , 


к 


and then the second equation gives 
кз (4 tan А — 3A)? = 4 (2tan A — ЗА) + 6 (4 tan À — ЗА) — 8 (1 + 4 tan A?) tan А 


= — 30А + 24 tan А — 32 tan АЗ, 
that is, 
32 tan АЗ — 24 tan A + 30A 
(4tan A — 3A)? 


Казр 


Hence, neglecting и, 


ә à 3— 15A 
cte УДАЙ р 13 16 tan ЛЗ — 12tan À + 


4tanÀ—3A 3 (4tan А— ЗА) 


It is extremely remarkable how little laborious this process is, as compared with any other 
likely to occur to a mathematician, for the elimination of e and v between the three equations 
3v--esinv PTS Lass ем v—esinv 


к= ———, n-= /[l—tanz, »v-AÀ-c4(-5—————. 
(1—ecos v)! 2/«I-e..2- (1—eoos v)! 


[3. The general expansion of к in terms of u] 


To calculate now the general expression for к. from the partial differential equation 


дк бк\ бк 
О= (-- эмо. 5+ 8tan А. 


We have* 


к= ук [0] 65; 


M usu i d wn 2 5 
= Eo [0] * – Хдл e LOD = Zo li- киа) [ОГИ = — оа; [01**; 
дк А бк уо (A 0-00-02 уе (A и 
= +) = (01 (А+) к; [0] (72 = Уе (А+ t) kig [0] "4 
= Удо (Ак; t ir) [0] 65; 
therefore, 
8 8 со , Ч oo 4.4 
+ 2{ап Асу + 35, = Xi [0] 6 (— к; + 2tan А (к! — кал) + ЗАки + 3i) = Xo 0; [0]. 
* [[0]-^— 1/;! See Mathematical Papers, Vol. 1, Appendix, Note 4, р. 468.] 
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The partial differential equation gives 
(ioo, [0]7* 6) x (250b, [0*5] = — 8tan A; 
that is, 


abo = — 8%ал à; 


and when ?> 0, 
Eiio tiv be [0]-6—? [0]-* = 0, 


а formula which may also be thus written* 


а= килку, b;-(3A—2tan2) к; + 2к tan À+ (32 — 1) к;; 
ар= к= ко,  bo= — ko + 2r tan A + (3А — 2 tan 2) к; 
and, since we know by the circular values of к and a that ку = ЗА, x, = 1, ко= 3, ку —0, we have 
а= —2, by—4tanA, 
values which satisfy, accordingly, the condition аб, = — 8tan A. Again 
Q,—kg—k,—kg, бу=2ку+ 2к|%ап А+ (3A — 2tanA) к, = (ЗА — 2tan А) «a + 2; 
therefore the condition 


Here 


(a+b) = ab, +аубу= 0 


becomes 
0= 2(2tanA— ЗА) x, — 4+ 4x, tan À, 
that is, 
374 cpi 
Ка 4tanÀ ЗА’ 


as before. Again, 
b,— (ЗА — 2tanA) кз --2k«;tan ÀA--bkaq, @=к.—ко, 


therefore the condition 
(a -- b); — a5 b, + 2a, 5, c agb, = 0 


becomes . 
0 — 4tan А («s — кз) + 2x3 (ЗА — 2 tan А) + 4к, + 2 (2tan À — ЗА) кз — 4«; tan А— 101, 
that is, | 
0= (4tan А — ЗА) кз — 4r tan А+ кз (ЗА — 2 tan А) — Зк,, 
in which 
016,8 ‚= 2(1+4%ап Аз) 
Ка Stani 3А’ '? (а%апА– 3А)? ` 
Therefore 
0 = 1к, (4%апА — 3A? + 4 tan А (1 + 4tan А?) + 2 (3A — 2tan Л) — 3 (4tan A — ЗА), 
and so | 
_ — 2 (16 ап А — 12tan À+ 152) 
o it (4tanÀ— ЗА)8 
as before. 


* Wo TA 
| @+у=а%+а + Ча +. 
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It might have been convenient to have calculated a,, b, as functions of A, and аз, 02 88 

functions of А and кз, before we proceeded to calculate кз from the equation 0 = (4 pb). ^? 
this way we should have had 

a= —2, bQ—4íanJX; 

2 4tanA 
hon diei Min апл SV 
j 2 (1+ 4tan А?) 
73 7 Ka — Ka 7 Кз (31an A ЗА)’ 
b, = (ЗА — 2 tan А) кз + 2x5 tan А + 5g 
16 tan АЗ — 36 tan À + 30A. 
= (ЗА — 2 tan A) кз — (&tanÀ— 3А) { 
so that the condition 0 = (a + b), — a4b, + 2a, b, + a,b, would have become 
1 
= -+ 2tan) 2) +4 8tanà3— 18% 15А 
0= 2к. {4 tan À — ЗА) + ББА 3i tan A (1 + 4tan А?) + 4 tan À + 8 tan À апл + 154}, 
that is, 
_ —2 {16%ап А — 12$ап À+ 15А} 
i Ti (4tan А— 3A)? , 

as before. 


In like manner, substituting this value of кз in аз and b,, we find 
la, (4tanA— ЗА)? = — (16 tan АЗ — 12tan À + 15А) + (&tan А2 + 1) (4tan A — ЗА) 
= 16tanA- 12A tan А? — 18А, 


therefore ; 
_ 4 (8$ап А — 6Atan А? — 9А). 
Ж, (4tan А— 3A)? + 
15, (4tan А — ЗА) = (2tan A — ЗА) (16 tan АЗ — 12tan А + 15А) 
— (4tan A — ЗА) (8 tan АЗ — 18 tan A + 15А) 
= 48 tan А? — 24A tan АЗ — 48A tan A, 
therefore 
_ 48tanA (2tanA — Atan à? — 2A) 
дб (4 tan А— 3A? 
Непсе 
ва 2325 _ 48 ап А (20tan А — 12A tan А? — 21А) 
ри (4tan À— ЗА)* Г 
Also 
аз + agb, = 4 tan À (e, — к) + 2 (2tan А— ЗА) к. — 4, tan À — 16s 
= 2 (4 tan A — ЗА) к, — 8 (кз tan А + 2кз); 
and 
1 (4 tanÀ— ЗА) к; = (4tan А+ 1) (16tan А? — 12tan А+ 15А) 
— (4 tan À — ЗА) (16 tan А* + 12 tan X? + 1) 
= — 80 tan АЗ — 16 tan À + 48A tan M + 96A tan А? + 18A, 
therefore 


‚ 12 (— 40tan 3 — 8tan А+ 24A tan А“ + 48А tan À? + 9А) 
T" (4tanA — ЗА) 
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Substituting these values in the condition 

(а- b), = ao b; + 3a, b; + 3a, b, + a3, — 0, 
we get 
— 1 (4tan À — 3А)5 x, = 3 tan А (80 tan АЗ + 36 tan А — 48А tan А — 108A tan А? — 392) 


+ 2 (64 tan M — 48 tan A? — 48A tan ЛЗ + 96A tan А — 452?) 
= 368 tan M + 12 tan А? — 144A tan А5 — 420A tan АЗ + 75A tan А — 905, 


that is, finally, 


къ йк ша ag б р 2 2 4 2— 9042 
ка dan À ЗД 4 tan А2 (92 tan А? + 3) + ЗА tan А (48 tan M + 140 tan А? — 25) + 9023]. | 


Thus if we only neglect (9, or (v — А), we have 


e З 16tanA?— 12tan А+ 15А 
4tanÀ—-3À 3' (4tan A — 3A)? 
14 — 4tan А? (92 tan А2 + 3) + 3à tan А (48 tan А + 140 tan А? — 25) + 902? 
8' (4tan А — ЗА)5 : 


к=ЗА-+ + 


[4. The results expressed in polar coordinates. ] 
Here, to recapitulate, А is the half of the angle 0, — 0, described in the time $, reduced to the 
case of equal final and initial radii vectores 7;, rọ; so that 
(£, — Eo)? + (n — то)? + (Ej — Eo)? = (7 + то) sin А?, 
(£j, эң, б, being the three rectangular heliocentric coordinates of a planet at the time t, and 
£o; 70, Со being their values at the time 0), and 


0,— 0\2 — (rj, —rgM? 0, — 09? 
р Т TE er эш. "о ыы у Yet 
sin À = (sin 3 | «(a 4d (os 2 | : 


¿is the following quantity, which is small for nearly circular orbits (such as here considered), 


ERE И: И 
té TEN 


ш being а constant such that the differential equations of relative motion are 
t= рат, qi ris G= eir’; 
ШР 
Vp (+ ту) 
8 being the principal function of the motion, or the definite integral 
t ( £2 /2 
sii | {= «Fat 
0 2 Ti 
which, when considered as depending on £,, т, £j, &0, 70, б and t, has (by my general theorem) 
its variation equal to 
£86, — 696, nio — mono + 080—8 +y- 80, 
2a being the axis major of the ellipse. 


and 
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А п")? E 
Accurately cos 2А = cos (0, — 0,) — 2 (л =) (сов а d 
therefore 
Р К —ту\? 0,— 04? 
sin (A — 3 (0,— 0,)} sin {А + 4 (6, — 69)) = (=) (cos | : 
Therefore, nearly, А = d more nearly 
0,70, 1 пм 0,— б. 
m1 aos mn 


more nearly still 


Ed. 4 Е] (cos М) cosec fo — 6+ : (=) cotan м 


2 +7 2\%+% 
" 39 + Е (; ^ 2 eoten ^ г % | Е 5 (=) opui 2 cotan (0,— в, 
and since, in general, 
(cotan x)? cotan 2x = аа старе = } (cotan 2? — cotan x), 


—т„\® 
we have, neglecting (i 4 
tt To 


_ 0,—0, 1 (r, — rg? 06,—0, 1 RR 0, — 0, ( 0—0)? 
А= ә sce) cotan 4 e (178 cotan 2 cotan № 


If we neglect the square of the eccentricity, then 


_ 0—0 E 2u 
Mea. к=0,— 0, - t | rj ny" 


Hb 


Tod. Tg. 


а= (8, AE 6,) м > To) + 


229 


And if, in this expression, we change r, into ry +7; — 7о, and still neglect the square of the ессеп- 


tricity, we find more simply (but less symmetrically) 
з = (0, — 60) V ur E; 
0 


because the coefficient of r, — 7, is equal to 


0, — bo к M (в cid Ey e 
а Мт, ar V. ° An) N 16 


and 6,— 0, —t + - is small of the order of the eccentricity. 
0 
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We can carry s to а further approximation as follows. We have 
M а= "ife. JT IN, 
accurately; therefore, nearly, 
8 Ti— To 1 (=) 
— =к41+ -z> t 
V ur | 4 2\4» 


neglecting the cube of the eccentricity. Neglecting the same cube, 


A NE 
5% + $ (" o) cotan " ; 
2ro 2 


А= 


l 
MICE TETTE CTS Rs 


where we have put » — V u[rà; 


E 2 _ ni(. Tiro)? ті 5 n) s s) | 
Ed м |26, А+ (1+ x ET ur ow iq ulmo 


therefore 


1 Злі ur 
2 10,—0, —nt 4- —- 
ме. "ы... п") , (1— ? (15м, dius n al йя: | жу]. 
Ж. Жыз АШ 4 | 2, а 16 * Stan Ini — Sn : 
and hence 
н NL APT 0 nt) + (^ TA + cotan ^H i ais 
Vas * e ir us 8 2 ` 8tan int—3nt 


4tan i (=. s (8, — 05 — nt) + 4 (6, — 6 — nt)? 
8tan E — 3ni j 


[5. The principal function in the three-body problem. ] 


Also, if we consider the system of the sun, M, and 2 planets, m, m, , the Principal Function S 
ofthe relative motion of this system mustsatisfy the two partial differentialequations following:* 


ss (tx) e) m) * Gc) 0а) c) Gn c.) 


х (5835 tss Б + ӨГ = mM ee тт, 
M lie дё. mòna 80,36, у + @ +@ Mon Ин é)? + (73 т)? + (£,,— Gy? ШР 


aC] э-нын) уе 


TM у + с mM. т,М mm 
М (8,8, бт бью 86086,0 


/ 


М3 + @ Е B TER, t 2, MU "7" éo)? + + (0,9— 79)" * (С, — £93) 6*5) 


* [These equations follow easily from (D*.), p. 153.] 
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If we put 
&=17,с080,,  7,—"rjSsinÓ, £,j—rgcos0,, «y—rosin6,, 
and similarly 
£,—7,0c080,, m,-r7,8in0,, E o=, Sino, m,o—r,osin6,g, 
we have 


55 9S в110,55 55 . „55 + 9088,55. 


artum, Mt de de" n 90) 
and therefore 


9,82 /8,8\2 [65\%_ 1/55\? eS? 6S? 1/59\? 11/59\? 

(sc) +5) -(&.) *alis) (iu + ) -(s.. +8 (у) , 
SAS 25 58 eS 8S 1 ôS ôg 159959 1 85868 
utente Herc giga.) 02 60+ (аи A * 


The first of the above two partial differential equations gives 


у ET E (891 
ôt 2Mm |Vàr, is 9C, 2Mm Ver, r 490, 9b, МУ +@ 


Mm, B as 85 | Е (RS INS pci 8S и) 
EZ n2 "n M — V)Orór,  rj90,r (90, M ór,r,80, ór,r,,90, 
m 65 65 тт 


о = , 


ТИШ Vid eris rr u08 (0,48) a 2° 


and the second may be similarly transformed. The first five terms may be made to vanish by 
employing elliptic values; the remaining terms give the perturbations. 

We shall assume that the inclinations are neglected, and shall put ,= = 0, {== 0. 
It is easy to perceive that in a system attracting according to Newton's law all the linear 
coordinates may be multiplied by any one common factor / (besides altering all the positions 
by any common motion of rotation), and all the masses may be multiplied by any other common 
factor l,, provided that the time # is multiplied by the factor ЙТ. And then, in the general 
expression for the principal function of relative motion, 


t t 
8= | т,й+ | U dt, 
0 0 


the coefficient U will be multiplied by l/l and T', will be multiplied by 1, ?/I3]71 = 12/1, and there- 
fore U dt and Т, й and finally S itself will each be multiplied by ЙЙ. Hence, in the present 
system, in which the orbits are both in one plane, 

8 Miri x fanct. (i Ab T, T, О.в, 0—60, 61—60 2 

rigorously, от, if we асое to express it so, then 


(+76 Ды. TAE Gaiti mm, 
"T ке. TE gs ОЛ Ит» Ун, 


W — V/r, туЁапс&. h, NN RE ce ию О, в), 


Tj Tg! Ta T, 
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к, А, ‹ having the same significance as before for the first orbit and к, , А,, г, similar quantities 
for the second orbit. In rigour, W involves т/М, т, [| M, but these small ratios may be treated 
as equal to zero in W, if we wish only to deduce those perturbations which are of the first order 
with respect to the disturbing masses. For greater symmetry we might put 


W= ja T Porn, 0,—0,-0,,— bo ЛА, 2. 1, WE ou ти? . 
2 2 Tt То Tatt, 


Thus, if we put for abridgment, as before, х = + (7, +7), и = М +m, and similarly 
a, = 5 (7,7,0), и, = М +т,; 


and if we put also 
Bus | Инт 1% 9 — Oc 8,o _ a+ 8 


ntr O Tatro 2 go 


we may write 


Е nu 


and may form two partial differential equations relatively to the function у as follows. 
S is explicitly a function of «, «, , À, À, , B, B, , « ‹, and 9, involving also the masses; its variation 
may therefore be thus expressed 


"D tmm, TE 4 (9, А, А,, ou. В, B, 


58 ,98 oS oS 85 8 
88 == Ba. ба ‚ + $у 9А+$ 5А, oM ^g 58+, ag, 4.58 E à ‚+$ 99 
Now 
— өл? 
sin А = (sin от "rx ej + (eos sm), 
and so 
. 1— f?sin (0,— 0o) 1 + cos (0, — 0,) 
Berry voor АЕ) Бына ПК 
Also 
Ла 2 cos? . 0,—0,)? віп А – В 
ир eet ~ уні Мает ч 
that is, 
1-+ сов (0,— 0%) cotanA  (1— B?)sin(0,—0,) _ ic At Eu 
sinJQÀ 1-88’  sin2 tok dendo 
therefore 
8A-1V I= Bë cosac Y (80-80) + cotan ру. 
7 
Now В= ЖП a= $ (rut ro), 
and hence 


 ón—óón В  1-f 1+ В 
В=—5 —– 5 (8+ дт) = p» 9—5 9. 


Also = “Ах 
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therefore E 
PIE m 1 |, 930-9 
И И cU fa С да. 
Непсе 
25 а 


dem p 858 deca e rette 
yd А+ A AE * s Ea Mj f 


&S 55\) (1—– B)8r,— (1+ B)àr, I- Bios (S. 90, — 96, 90, 
s iam) EIS prem" 
that is, 
ДМ Ды (б-ны ы, шч санасы, 
al Е" "ET 0x 2 « à 2 би 2 а, ё 2 


ГА 


85 , Beotana 858 а 3 BiootenA, (5 - 5] 
Ый т шү м: ТОЕ Г И V. S 
(15 Ам шй ык ВЕ coe (095. PEE 905 Ут Bi cosec à 
x a, --V/1— В? cosec А a IN В? cosec à? 
x ый 50—50» 85 50,,—80,-- 80, — 80, 

hk) А ТЕ 


From this equation we immediately deduce the partial derivatives of S with respect to the 
variables f, т, т, 6,, Ө, ro, 7,0, 09, and 0,9, and we have then to substitute these values in the 
two partial differential equations (for the case of null inclinations) 


_ 88 no os) Gn M 4m, (( 8S? 85 у-и 
© 2Мт Ver, 2) |+ 2Mm, (2) * (as т; Tu 


cos (8,,— 6) coe 88 va sin (0,,— 6) ts 8S 59 88 ) 


4 


M бт, бт ШЕЛ M r,r,90, ӛт, т, 100, 
тт, 
Vh. ri — 2r ricos (0, — 6) 
and 
о= + эии s) + 55 \ Mtm, x) 58 тн 
EA ôt 2Mm ёт AS B 2Mm, ls. Hur; То 7,0 
UA (S S 595 68 juo (E 58 в 
M ror, тобот, 80, M Oro одб Этот, 950, 
mm, 
125 +72 — 2r, org сов (8,9 — 60) 
We have : 
0,—0,, 0,—0 0,—0,, , = 
M DR ub ее 0,707920. 7, 
and 
0,— 66 cosà . 6—6, sin À MEI 
сов——ь УЕ sin hi m 1 — f? cosec 
0,—0,9 cosà . 0,—0,9 us VET oec X: 
cos —— 9 = , mi dui Де” ! V 1 — В? совес A7; 
dE dp MD v ME 


HMPII 30 
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therefore 
cos (0,,—6,)) = К cos 9 — Lsin9, вір (0,,—0,) = К cos + Lsin9, 


where we have put for abridgment 


E- {ча E ME 6 ®) сов А cos А, + sin Asin А, V1 — В? созес X2 V 1 — В, совес Xj 
сов [74,8 — 179 | С=с, E у Р. STEEP 


2 V1—Bv1- В? 
Liin (муе о &— ^)- cos Asin A, V1 — SET x EN М/1 — В? cosec X? 
| 2 2437 ул рут р у 


Thus, making the above substitutions, the first partial differential equation becomes rigorously 
(for the case of null inclinations) 


ш8 1 Mm Mm, 


=} [498,1 [n88. Mm . Mm, 

ee DENT «(1+8) а, (1+4,) 
ааа да ААА За о А 
Е 90а, (1+ B) (1+ B,) (К cos 9 — L sin 9) 


, 


+з |$ T 2 204 Lm Ee (m x) 
ыа“, КВ, #0, TR (5, в] 
METETE e (0а све) ва 

4 a sui» 1 — B? cosec А? э, -5.) sl 
editt „а-а-а иШ М) 


[ева +1-В$я-- ы, heri в ш) 


xui Lsin9 [роже (88-18) 88] 


"Мәә (1+ B) (1 4 B;) 5 &) 89 
——MX———,í(98 8S8 5S) Ksin9-4.Lcos9 (|, [98 88, 858 
di P Mid -&.)*3] "AMax, (+В) (13 B) [Vic Brosse (SS 5) - | 
55 , B, cotanà, 58 88 
d gp BJgg - А, ++ А, (5. -&)| 
Ait cent) 
‚88 58 + Beotana 8S 8 
х |а +а- Bg - TP: Gi- T l 


and the second partial differential equation may be rigorously formed from this by merely 
changing the signs of В, B, , 9. 
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[6. Approximations in the case of nearly circular orbits.| 
If we entirely neglect B and В, , or rather suppose them rigorously to vanish, then the partial 
differential equations become 


o 1 [№85 M Мт мт. ы „тт. 
ELIT oj ё; « а, Мо + о — аа, cos (9 -- 4, — À) 
М+т | 58 88 ‚ М+т, | 88 | 98 8) 
BEAN &.*ap 10 0% | +} | а +8, В+) т 
EE М+т, (88 88 às, 
*&MmoRV8A & 55) ^ BMm, о (5А pae 
пова) at $8 +58 oS 85 as 
taa а -i0- "EHI ба В. # (А, +e r3 
cos (8+, —A) [85 8S BS, (88 8BS,58 
ах Ае WV S A 39 
sin (9-2, made ota 88 28 
ху IM M P tag i^ + | 
 Sin(9--A,— —А) (55 98S 85 ТА 85 85 
VET M оа (5 5, 2 tag +) zi 


together with the other similar partial differential equation derived from this by changing the 


55 
sign of 9 and those of x в, '8B,' . In these equations no power of any mass is neglected nor any 


of г, 1,5 but if we put NE 


® «, и“ 
S Mme |+ Mn, v, | + тт, : 
ш №, PAM 


and neglect тт, and mm? , the equations become, after being reduced by the partial differential 
equations of elliptic motion and divided by mm, /16, 


о- [А [ir («- AS- (з 2x aF- ^d 

die ee etenim nt 
+ NE С-КЕ" ИР үт B- (+) 
e Fez t ET 


8 д дк, дк, 
5n ms E a-a g) (к, M ENT 2d 


s Spa. 
ШООК T4 дк дк бк, дк, 
(ams) (в) 
2M бк, бк, м) -( ^ бк _ x) бк, x, 
TUER віп(+ 9+2, — E Nis ide EI 5, к 2x 3, x А S 


16 


"ON E CM ПЕНИИ НИНЫ i^ MU m а 
Ма? + a3 — Зоо, cos ( + 9+ А, —А) 
30-2 
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This double equation is rigorous so far as powers of « are concerned, though В, В, and тт, 


т?т, have been neglected. But if we change к and к, to their approximate values ЗА + and 
ЗА, + ‹,, and then neglect г, г, , the double equation becomes 


MER p, nA + A- je, [х+«, а Gs Ld Mtr 
72x М \$А 59 E 8A, 89 М ua V px, 
1 алы Eu mI 
Маз od — Эа, cos ( + 9 +А, – А) 
If we now put 
| В, 
and therefore А 
— —- A A, M XX 
pen 2 M 2 — mm — e a " 
pax = na Bx, —nm, = я,’ Vu Vx 
and make as a sufficient approximation и/М = p, | M = 1, we shall change the double equation 
to the following: 


ом А (Se) es (20 = 


2 
that is, 


о, 
У аА ЁЛ, сов ( + 94A, —А) 
o лы 


Y БА ze) cos (А, — À) cos 9 — ——— 


REN LLLA cen 
ГА E 
«y nen 


, 


+ ЭХЕ cos ( + 9--À, —5) 


1 
ACUERDO 
ЛЕЙ; С ы à,-X 
is (i) — А; сов (— 9+ , 
апа { 
a ү 5 AM, 1 
ve dd E E sin (А, — A) sin 9 — 


m oa BEBMINIC Nt ome erga ОРДЕ 
А «(D ое +À, -2) 
JG) АЖ УА. 


1 
—————7* 
-4 л\- w 
hh) -3(z) eos (— 9 4- A, А) 
$ M T 
We have seen that =% =[®) ; if then we put for a moment a/*, 7 * and А, –А= А, we shall 
have : 
А : x a" № 
— = 4 = — йш —— у 
À & ЕЛ А a3 -—1' + «ё 1 


MOM 204XÀ ар M 2500, 
5. &' 35 M^ DÀ 80 ELA А’ 
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and hence 
ор. AW, 1_А бу 
ÖA, А,дА А, SAY 
Therefore 
MS an à 2 1 1 
02 V1-4« (1—a'3) — + —  cosA'cos 9 — —— _———== rr RM romane cU 
prse Py Jirar acos А) М1 ж+—да'оов(—9+Х) 
and 
— Ww x 1 1 
0=V1 +a (1—a'3) — — —— вір Авію 9 – аа а 
( 59 Уо М1 -+«`?— 2a cos (9--X) Уна Эа‘ со(—9+Л\) 
Treating «' as constant, this gives 
е Зы ‚ 89 + 8A 
М + « (1 а) 8 2 — (sin A sin 989 — cos А“ cos 9.81") + = 
di Va ар М1 +072 — 90 cos (9 А) 
5А` — 89 


+ аа Я ? 
V/1-- 0% — 20 cos ( 9 -À') 
that is, 


MEAS 94^ ` 
V 1-4. (1—0) = funct. («`) -Asin À` cos 9 «[ Md — 
c 


9—» М1 +2 — 2x' cos TT 


But the arbitrary function of «`, introduced by this integration, must be identically zero because 
otherwise it would not vanish with à`, as it must do, since y does so, independently of «` and 9. 
Thus, finally 


SCR lac. кА nad 1,$ 944-24 dó" 
Su vs VE T E - 3 pe 24 в 
ata, Jj — Va, М1 + a y ` ES cos 9 sin (А, — А) + io" on P Vas 3 o3 бах cos 8 
which may also be written 
(+, 4) 3 | | | CA d6" 
- -2josPsnQ,-34| áp me» 
TIE EM оов 9 sin (А, AE | аА ACE аА ФА oos 0) 


ог in the alternate form 


«T 


} 94 4 3 AM -i 
а КЕ — À) + NO +( ) EL dé. 


